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We present a detailed theoretical investigation on the magnetic properties of small single-layered 
Fe, Co and Ni clusters deposited on Ir(lll), Pt(lll) and Au(lll). For this a fully relativistic ab- 
initio scheme based on density functional theory has been used. We analyse the element, size and 
geometry specific variations of the atomic magnetic moments and their mutual exchange interactions 
as well as the magnetic anisotropy energy in these systems. Our results show that the atomic spin 
magnetic moments in the Fe and Co clusters decrease almost linearly with coordination on all three 
substrates, while the corresponding orbital magnetic moments appear to be much more sensitive 
to the local atomic environment. The isotropic exchange interaction among the cluster atoms is 
always very strong for Fe and Co exceeding the values for bulk bcc Fe and hep Co, whereas the 
anisotropic Dzyaloshinski-Moriya interaction is in general one or two orders of magnitude smaller 
when compared to the isotropic one. For the magnetic properties of Ni clusters the magnetic 
properties can show quite a different behaviour and we find in this case a strong tendency towards 
noncollinear magnetism. 



PACS numbers: 75.75.-c, 75.75.Lf, 75.70.Tj, 75.70.Ak 

I. INTRODUCTION 

The magnetism of surface supported clusters has been 
the subject of intense research activities over the last few 
years as such systems often show peculiar and unexpected 
magnetic behaviour. These exceptional magnetic proper- 
ties arise from the reduced dimensionality in combination 
with spin-orbit coupling (SOC) which can cause complex 
interactions among the atomic magnetic moments. In 
this context clusters of magnetic id transition metal el- 
ements deposited on 5d noble metal substrates are very 
interesting as for these systems spin-orbit driven effects 
mediated by substrate atoms with large SOC are most 
prominent. With technical or chemical applications in 
focus, there is a growing need to understand the trends 
and principles behind the manifold of magnetic properi- 
ties for different cluster and substrate materials as only 
this will make it possible to anticipate which magnetic 
properties may result from a particular cluster/substrate 
combination. 

In previous experimental and theoretical investigations 
on the magnetism of atomic clusters on surfaces it was al- 
ready demonstrated that their magnetic properties differ 
strongly from the magnetic properties of the correspond- 
ing bulk materials and that this has its main origin in the 
reduced atomic coordination of cluster sites which in fact 
has a strong impact on the local spin and orbital mag- 
netic moments^— . More recently it was also shown that 
for 3d clusters or monolayers on 5d metal surfaces SOC 
induced effects on the spin configurations also play an im- 
manent role causing various noncolliniear magnetic struc- 
tures^—. This SOC induced noncolliniear magnetism is, 
however, intrinsically different from the spin frustrations 



that may arise e.g. by a competition between ferro- and 
anti-ferromagnetism or that may be present in systems 
where the magnetic and geometric symmetries are incom- 
patible 7 -^. 

Unfortunately, each of the theoretical studies pub- 
lished so far were aimed at only one or two combinations 
of the cluster and substrate materials and often only very 
few cluster sizes and shapes were investigated. In addi- 
tion to that comes the fact that many theoretical inves- 
tigations have focused only on some selected magnetic 
properties as for instance the magnetic moments and ex- 
change interaction but leaving out important informa- 
tion concerning the magnetic anisotropy energy (MAE). 
Moreover, due to limitations which are present in all the- 
oretical schemes it is often also problematic to compare 
results obtained for different systems by different groups 
which use different methods. Thus, in order to obtain a 
more complete picture about the trends in the magnetic 
properties of deposited clusters one needs a sufficiently 
large self-contained set of results for interrelated systems 
which are obtained by the same method. This motivated 
us to calculate a large spectrum of the magnetic proper- 
ties for sets of Fe, Co, and Ni clusters of 1-7 atoms on 
Ir(lll), Pt(lll), and Au(lll) surfaces, within a unified 
fully-relativistic Green's function formalism. Moreover, 
we studied also complete monolayers as reference systems 
for the sequences with increasing cluster size. This en- 
ables us to analyse a large pool of data which are directly 
comparable because they were obtained by the same pro- 
cedure. We found that the magnetism of Fe and Co clus- 
ters on all investigated surfaces follows common patterns 
that can be understood by considering the coordination 
numbers of atoms in the clusters and the polarisability 
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of the substrate. For Ni clusters the situation is more 
complicated and some of the systematic trends observed 
for Fe and Co clusters are absent. 



II. COMPUTATIONAL FRAMEWORK 

The calculations for the investigated cluster and mono- 
layer systems were done within the framework of spin 
density functional theory (SDFT) using the local spin 
density approximation (LSDA) with the paramctrisation 
given by Vosko, Wilk and Nusair for the exchange and 
correlation potential^. The electronic structure has been 
determined in a fully relativistic way on the basis of the 
Dirac equation for spin-polarised potentials which was 
solved using the Korringa-Kohn-Rostoker (KKR) multi- 
ple scattering formalism 11 ^. The calculations for surface 
deposited clusters consist of two steps. First the host sur- 
face is calculated self-consistently with the tight-binding 
or screened version of the KKR method^ using layers of 
empty sites to represent the vacuum region. This step is 
then followed by treating the deposited clusters as a per- 
turbation to the clean surface with the Green's function 
for the new system being obtained by solving the corre- 
sponding Dyson equation 12 . This technique avoids the 
spurious interactions between clusters which may occur 
if a supercell approach is used instead^. 

For all systems discussed below the cluster atoms were 
assumed to occupy ideal lattice sites in the first vacuum 
layer and no effects of structure relaxation were included. 
The substrates were simulated by finite slabs which con- 
tained 37 atomic layers and we used lattice parameters 
of 3.84 A, 3.92 A and 4.08 A for Ir(lll), Pt(lll) and 
Au(lll), respectively. The surface calculations were con- 
verged with respect to fc-point integration. For the sur- 
face Brillouin zones a regular fc-mesh of 100 x 100 points 
was used which corresponds to 1717 fc-points in the ir- 
reducible part of the Brillouin zone. The effective po- 
tentials were treated within the atomic sphere approxi- 
mation (ASA). The occuring energy integrals were evalu- 
ated by contour integration on a semicircular path within 
the complex energy plane using a logarithmic mesh of 32 
points. The multipole expansion of the Green's func- 
tion was truncated at an angular momentum cutoff of 
dax = 2. For selected surface and cluster systems calcu- 
lations with £ max = 3 were also performed which showed 
that this causes a more-or-less uniform increase of the 
local spin moments by 3-5% and of the local orbital 
moments by 3-10%. This indicates that the systematic 
trends in the spin and orbital magnetic moments are well 
described by ^ m ax = 2. 

For the representation of the interatomic exchange in- 
teractions we made use of the rigid spin approximation 13 
and mapped the magnetic energy landscape E({ik}) onto 
an extended classical Heisenberg model for all atomic 
magnetic moment directions {e^}. The corresponding 



extended Heisenberg Hamiltonian has the for m" 15 : 

~2 J2 ®v ■ & x ^ +J2 K ^ ii "> ' w 

where, the exchange interaction tensor has been decom- 
posed into its conventional isotropic part Jy, its trace- 
less symmetric part J$ and its anti-symmetric part 
which is given in terms of the Dzyaloshinski-Moriya (DM) 
vector Dij. We calculated the coupling parameters 
and DM vectors following the scheme by Udvardi et 

The anisotropy constants Ki{e{) account for the on-site 
magnetic anisotropy energy associated with each individ- 
ual magnetic moment oriented along e,. The magnetic 
anisotropy energy AE is usually split into two parts, the 
SOC induced magnetocrystalline anisotropy AE soc and 
the so-called shape anisotropy AEdd caused by magnetic 
dipole-dipole interactions, i.e. 

AE = AE soc + AE dd . (2) 

AEdd can be determined classically by a lattice summa- 
tion over the magnetostatic energy contributions of the 
individual magnetic moments or in an ab-initio way by 
using a Breit Hamiltoniar>i£. Here, we used the classical 
approach to calculate A-Edd for the full monolayers while 
we found that for clusters containing just a few mag- 
netic atoms AEdd is negligible. The magnetocrystalline 
anisotropy energy Ai? soc was extracted from magnetic 
torque calculations which are described in more detail in 

As discussed recently by Sipr et alJ£ the approxima- 
tions and truncations mentioned in this Section result in 
a limited accuracy concerning in particular the values of 
AE soc . However, this does not hinder our analysis of the 
general trends of Ai5 soc with respect to cluster geome- 
tries as well as different cluster/substrate combinations. 

III. RESULTS AND DISCUSSION 

A. Magnetic moments 

Fig.[T]shows the considered cluster geometries together 
with the calculated values of local spin (/z sp i n ) and orbital 
(Morb) magnetic moments for Fe clusters of 1-7 atoms as 
well as full Fe monolayers deposited on Ir(lll), Pt(lll) 
and Au(lll). For identical Co and Ni clusters the cor- 
responding data are presented in Figs. [Hand O In addi- 
tion, these figures also show the induced spin magnetic 
moments of the respective substrate atoms that are ad- 
jacent to cluster atoms. One can see that in some cases 
there are considerable variations of [X S pm and /i or b be- 
tween the different sites of the deposited clusters. The 
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FIG. 1: Cluster geometries for Fe clusters of 1-7 atoms supported by Ir(lll), Pt(lll) and Au(lll). The local spin and orbital 
magnetic moments at Fe sites are given by the upper and lower numbers, respectively. The spin magnetic moments for nearest 
neighbour substrate sites are also shown. The data presented within frames give the corresponding monolayer values. 



magnetic moments depend not only on the position of 
the site with respect to other Fe, Co or Ni atoms but 
also on their position with respect to the underlying sub- 
strate atoms. This can be observed for example when 
inspecting the two differently located compact trimers or 
the cross-shaped five atom clusters in Figs. [TJ [5] and 
Clusters supported by Pt(lll) have largest /J, sp i n when 
compared with Ir(lll) and Au(lll) while the /^ OT b val- 
ues are increasing from Ir to Pt to Au. 

There is a big difference between the induced mag- 
netic moments in the Ir(lll) and Pt(lll) substrates on 
the one hand and the Au(lll) substrate on the other 
hand. Ir and Pt atoms which are nearest neighbours of 
any Fe or Co atom have a relatively large fj, sp i n of up 
to 0.15 fj,Bi while corresponding Au atoms have always 
small negative fj, sp im not larger than 0.03 /xb in the ab- 
solute value. Substrate atoms with a larger number of 
Fe, Co or Ni neighbours usually have a larger /j, sp i n than 
substrate atoms with a smaller number of neighbouring 
cluster atoms. However, this is not a general rule as seen 
for the substrate atoms adjacent to the central atom of 
the cross-shaped Fes and C05 or to the differently located 



compact Fe3 and C03 clusters on Ir(lll) and Pt(lll). 

The orbital magnetic moments induced in the sub- 
strate atoms are always small: they can reach up to 
0.03 /is for Fe and Co on Pt(lll) while being smaller 
than 0.007 /is for Ir(lll) and smaller than 0.004 /zb for 
Au(lll). Except for the Au(lll) substrate atoms /i or b 
is found to be always parallel with /i S pm- The finding 
that Pt is the most polarisable of the three elements and 
that Au is less polarisable than Ir is consistent with ear- 
lier theoretica l 20 i 21 and experimental works 2 ^— for mul- 
tilayer systems. This high spin polarisability of Pt can 
be ascribed to its high spin susceptibility that in turn 
is caused by its relatively large density of states at the 
Fermi level leading to a large Stoner product (see below) . 

By plotting the local magnetic moments as a func- 
tion of the coordination number one can visualise the 
site-dependence of /i S pin and fj, ib- Such plots are shown 
in the insets of Fig. [3J where only neighbouring cluster 
atoms are considered in defining the coordination num- 
ber. Sites with a lower coordination number generally 
have larger ^ sp ; n and /x or b than sites with a higher coor- 
dination number, with Ni on Ir(lll) being the only ex- 
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FIG. 2: Cluster geometries for Co clusters of 1-7 atoms supported by Ir(lll), Pt(lll) and Au(lll). The local spin and orbital 
magnetic moments at Co sites are given by the upper and lower numbers, respectively. The spin magnetic moments for nearest 
neighbour substrate sites are also shown. The data presented within frames give the corresponding monolayer values. 



ception. For Ni on Fr (111) the magnetic moments quasi- 
oscillate strongly with changing cluster size or shape and 
we find that /i sp i n for the adatom (0.38 /xb) is smaller 
than for a Ni atom in the full monolayer (0.49 For all 
other cluster/substrate systems considered in our study a 
quasi-linear relationship between [i sp i n and coordination 
number is found. Interestingly, increasing the coordina- 
tion number for the atoms of such small clusters leads to a 
stronger reduction of [i S pin when compared to equally co- 
ordinated atoms in larger clusters or full monolayers. For 
example the central atom of a compact 7-atom cluster has 
always lower /x S pin than a monolayer atom. One can also 
see that the corresponding orbital magnetic moments are 
much more sensitive with respect to coordination than 
the spin magnetic moments. While the insets in Fig. [4] 
show a strong decay of /i or b with increasing coordination 
for Fe and Co clusters on all three substrates the orbital 
magnetism in Ni clusters behaves non-monotonously. 

An analysis of the average spin and orbital magnetic 
moments as function of cluster size is shown in Fig. 2] 
For the three and five atom clusters the lower /i sp i n and 
Morb values correspond to the compact clusters. All clus- 



ters have largest ^ sp in when deposited on Pt(lll) fol- 
lowed by the Au(lll) substrate. The lowest /u sp in values 
are obtained for deposition on Ir(lll). The highest val- 
ues of /iorb) however, are found for clusters deposited on 
Au(lll) where the interaction between cluster and sub- 
strate atoms is weak and the lattice constant largest. 

Concerning the trend of /i sp i n for the three different 
substrates there are two competing effects that must be 
considered. At first there is an increase in the lattice con- 
stant when going from Ir (a=3.84 A) to Pt (a=3.92 A) to 
Au (a=4.08 A), i.e. as the atoms of the clusters occupy 
ideal lattice sites, their distance from the substrate is 
largest in the case of Au. This means also that among the 
discussed substrates, the interaction between adatoms 
and substrate is smallest for Au and one would therefore 
expect that clusters deposited on Au(lll) would have 
the largest /Lt S pin values. On the other hand hybridisation 
of the electronic states between adatoms and substrate 
leads to a small charge transfer of minority 3c?-electrons 
from the cluster atoms into empty 5(i-states of adjacent 
substrate atoms, thereby increasing fj. sp i n for the clusters. 
This however, happens only for the spatially extended 
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FIG. 3: Cluster geometries for Ni clusters of 1-7 atoms supported by Ir(lll), Pt(lll) and Au(lll). The local spin and orbital 
magnetic moments at Ni sites are given by the upper and lower numbers, respectively. The spin magnetic moments for nearest 
neighbour substrate sites are also shown. The data presented within frames give the corresponding monolayer values. 



5<i-states of Ir and Pt with their 5d-states having an ap- 
preciable energetic overlap with the minority 3d-states of 
cluster atoms. This can be clearly seen from the den- 
sity of states curves which are presented in Fig. [5] In 
contrast to this there is almost no or only little interac- 
tion between the minority 3d-states of cluster atoms with 
the energetically low-lying 5d-states of Au. Besides, the 
hybridization between the cluster-derived and substrate- 
derived states leads to energy lowering of the Fe, Co and 
Ni 4p-states and this lowering is again more pronounced 
for the Ir and Pt substrates than for the Au substrate. 

This causes an additional charge redistribution within 
the cluster atoms, i.e. 4p-states become occupied, again 
at the cost of minority 3(i-states. In this way Fei de- 
posited on Pt(lll) ends up with about 0.1 electrons less 
in the minority 3c?-orbitals and thus a slightly larger spin 
magnetic moment when compared to the deposition on 
Au(lll), with this effect becoming stronger for Co and 
Ni. Sipr et al.— have demonstrated this behaviour for 
Co clusters on Pt(lll) and Au(lll) by performing cal- 
culations for Co clusters on an Au(lll) substrate using 
the lattice constant of Pt. This also showed that the ob- 



served increase in /i or b can be solely attributed to the 
larger lattice constant of Au. 

B. Density of states 

In Fig. [5] we show the spin-resolved density of states 
(DOS) for the adatoms, dimers and the central atoms 
of the 7-atom clusters as well as the corresponding full 
monolayers. The DOS of the respective topmost atomic 
layer of the clean substrates is also shown by the brown 
areas together with the respective Bloch spectral function 
A B (k, E) along the high symmetry lines r-K-M-f in the 
two-dimensional Brillouin zone presented in the first col- 
umn. A B (k,E) can be interpreted as a fc-resolved DOS 
revealing the detailed features in the electronic structure 
of the three different substrates. The large grey-shaded 
regions in the A B (k,E) diagrams in the first column of 
Fig. [5] correspond to electronic bulk states of the under- 
lying substrates, while the sharp black lines represent 
surface states localised within the topmost atomic layer 
of the clean surfaces. The blue and red regions arise from 
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hybridisations between surface and bulk states. Clearly 
visible are for instance the Rashba split surface states 
around f for Pt(lll) and Au(lll). 

For Ir and Pt there is an appreciable energetic overlap 
between electronic states located at the substrate and 
states located at cluster sites resulting in hybridisation 
with a prominent broadening in the cluster DOS. The 
situation is different for the Au substrate, where the en- 
ergetically low-lying states of Au can only hybridise with 
the majority states of Fe while for minority states of Fe 
as well as for all states of Co and Ni, there are no en- 
ergetically close Au states to hybridize with and hence 
very distinct atomic-like features prevail in the DOS of 
cluster atoms in that case. With increasing number of 
cluster atoms a complex fine structure appears in the 
DOS which also broadens appreciably with increasing co- 
ordination numbers of the cluster atoms. Moreover, the 



presented DOS curves for the central atom of the 7-atom 
cluster demonstrate that the DOS of deposited clusters 
acquires very quickly the main features which are present 
in the DOS of a complete monolayer. However, as for 
such small clusters the DOS at the Fermi level varies 
strongly with changing the number of atoms so do the 
corresponding chemical and magnetic properties in this 
finite size regime. 

The decreasing overlap between states located in the 
substrate and located in the clusters when going from Ir 
to Au explains the finding that /i or b is largest for clusters 
on Au where the atomic-like character of the DOS pre- 
vails. In this context, however, the size of fi sp i n cannot 
always be directly related to the overlap between cluster 
and substrate DOS as this overlap is smaller for Au(lll) 
than for Pt(lll) and yet /%,in is largest for clusters on 
the Pt(lll) substrate. 
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In the same way also the induced magnetisation within 
the substrate depends on this mutual energetic overlap 
between 3d and 5d-states, which explains the very small 
induced magnetic moments in the case of Au(lll). But 
also here, one should keep in mind that the polarizability 
of the substrate atoms is determined by the Stoner prod- 
uct / • Np with the exchange integral / and the number 
of states at the Fermi level Np (I-N F = 0.29 for Ir, 0.59 
for Pt and 0.05 for Au)2£. 



C. Exchange Coupling 



The calculated isotropic exchange coupling constants 
Jij are presented for all clusters (except for the cross- 
shaped 5-atom and 6-atom ones) in Tab. |TJ Positive and 
negative values of J,j correspond to ferromagnetic and 
antiferromagnetic coupling, respectively. In Tab. [TT] we 



show the sum of all couplings related to a particular atom 
at site i. 



(3) 



The effective exchange field J 1 can be seen as the total 
strength by which the magnetic moment at site i is held 
along its direction by all other atoms. 

From the data given in Tab. Q] one can see that Fe and 
Co clusters show a strong ferromagnetic nearest neigh- 
bour coupling in the range of 40-140 meV while for Ni 
clusters the Jy values are often much smaller. In the case 
of Fe and Co the couplings between nearest neighbour- 
ing atoms are about one order of magnitude larger than 
couplings between more distant atoms, i.e. the coupling 
strength falls off very rapidly with increasing the inter- 
atomic distance. For Ni clusters, however, and especially 
for Ni on Ir(lll) where the couplings are very weak this 
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TABLE I: Isotropic exchange coupling constants Jij in meV 
for Fe, Co and Ni clusters deposited on Ir(lll), Pt(lll) and 
Au(lll). The icons in the left column indicate the corre- 
sponding cluster geometry as well as the cluster sites i and j, 
respectively. The last line gives Jij for nearest neighbouring 
sites within the full monolayer (ml). 
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trend is less pronounced. These results also show that 
there is an occasional weak anti-ferromagnetic coupling 
between more distant atoms for Fe clusters, which how- 
ever, gives only an insignificant contribution to the cum- 
mulative J 1 of each respective atom. As each Jij contains 
by definition (see the Heisenberg Hamiltonian in Eq. (Q}) 
the product between the involved spin magnetic moments 
/is pin and /ig pin the coupling is, naturally, largest for Fe 
and smallest for Ni clusters. Moreover, the nearest neigh- 
bour exchange coupling among the Fe, Co and Ni cluster 
atoms is larger than our corresponding values for stan- 
dard bcc Fe (37.8 meV), hep Co (26.3 meV) and fee Ni 
(4.8 meV). 

Apart from the magnitude of the spin magnetic mo- 



TABLE II: Effective exchange field ,F in meV for Fe, Co and 
Ni clusters deposited on Ir(lll), Pt(lll) and Au(lll). The 
icons in the left column indicate the corresponding cluster 
geometry as well as the cluster site i. The last line gives J* 
for the full monolayer (ml). 
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ments also atomic coordination as well as substrate ef- 
fects play an important role. Especially for Fe clusters 
the values between low coordinated cluster atoms are 
often much larger when compared to atoms with higher 
coordination. Nevertheless, the effective exchange field 
J 1 increases monotonically with increasing coordination, 
i.e. given a fixed number of Fe, Co or Ni atoms the most 
compact structure will form the most stable ferromag- 
netic state. 

As one can clearly see from the data in Table Q] in com- 
bination with the cluster geometries given in Figs. [T][3] 
the isotropic exchange coupling is also affected by the 
arrangement of cluster atoms with respect to the under- 
lying surface sites. Looking at the two different compact 
Fe and Co trimers on Ir(lll) and Pt(lll) for instance we 
find that the coupling values differ by about 8-10%. For 
the 7-atom Fe cluster on Ir(lll) and Pt(lll), however, 
for nearest neighbouring edge atoms varies by as much 
as 20-45%, respectively, whereas the corresponding cou- 
plings for clusters on Au(lll) do in general not exhibit 
such a pronounced dependence on the atomic position 
with respect to the substrate atoms. The latter seems 
also to be the case for Fe and Co 7-atom clusters with an 
identical configuration on Cu(lll). This was studied in 
detail by Mavropolus et al<£ and the Cu substrate atoms 
also do not seem to participate in the exchange coupling 
of the Fe and Co cluster atoms. Therefore, we ascribe this 
substrate effect to the large spin-polarisation within the 
the Ir and Pt surface atoms while the weak induced mag- 
netism in Cu and Au causes only minor variations in the 
exchange coupling of equidistant cluster atoms. These 
irregularities in the couplings underline that transferring 
Jij coupling constants obtained from bulk calculations to 
low-dimensional finite nanostructures will lead in general 
to unreliable results. 

For Fe and Co clusters the magnitude of the isotropic 
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exchange interaction is quite similar for all three inves- 
tigated substrates. Ni clusters, on the other hand, have 
comparable nearest neighbour J^- values only when being 
deposited on Pt(lll) and Au(lll) while deposition on 
Ir(lll) reduces the coupling strength to just a few meV. 
This results in a quite small effective exchange field J l per 
Ni atom in the order of about 10 meV. As the exchange 
interaction is so small for these Ni clusters, there is a 
pronounced tendency that their magnetic ground state 
deviates strongly from a collinear configuration (see be- 
low). 

The coupling of magnetic cluster atoms to the induced 
magnetic moments in the substrate is always very small. 
Jij is about 2 meV between Fe or Co cluster atoms and 
topmost layer atoms of an Ir or Pt surface. The small 
induced moments in the Au(lll) substrate couple anti- 
ferromagnetically to the cluster atoms. Here the nearest 
neighbour J^'s are only in the order of 0.1 meV being 
of similar magnitude as the ferromagnetic coupling of Ni 
cluster atoms to Ir or Pt surface sites. 

In addition to the isotropic Jy coupling con- 
stants Tab. IIIII shows the complementary data for the 
anisotropic exchange interaction. For clarity we present 
here only the magnitude of the DM vector \Dij\ which 
can be seen as a measure of the driving force towards a 
non-collinear spin configuration. Given the fact that the 
SOC strength is comparable in Ir, Pt and Au one can see 
from the data in Tab. IIIII that for any given cluster there 
are often strong variations (without any clear trends) in 
\Dij\ upon deposition onto different substrates. As dis- 
cussed above for the Jy values, we find here an even more 
pronounced dependence of |-Dy| on the position of clus- 
ter atoms with respect to the underlying substrate atoms 
and the results show in addition that the relative decay 
of the DM interaction with increasing interatomic dis- 
tance is much weaker when compared to the corresonding 
isotropic exchange coupling. 

Albeit that \Dtj \ is between one or two orders of mag- 
nitude smaller than the isotropic exchange coupling, it is 
not negligible. For Fe n on Ir(lll) we obtain a relatively 
strong DM interaction which is in accordance with the 
recent findings of Hcinze et al£ and of von Bergmann 
et ali^Z, as well as Deak et al^ for Fe/Ir(001) which all 
demonstrate that these systems show a strong tendency 
towards non-collinear magnetism. Moreover, our results 
also show large \Dij | values for Fe„ and Co„ clusters from 
which we conjecture that this may also lead to complex 
magnetic structures within extended Fe and Co nanos- 
tructures on these substrates. In fact the sometimes 
experimentally observed, unexpected low magnetic mo- 
ments in Fe- and Co-Pt(lll) systems may be caused by 
this mechanis m 29 ' 30 . For Ni clusters the DM interaction 
is always very important with respect to the isotropic ex- 
change coupling as both quantities are often of the same 
order of magnitude. Thus, one can expect the presence 
of non-collinear magnetic ordering in Ni clusters on all 
three substrates. It should be stressed that this non- 



TABLE III: Anisotropic exchange coupling parameter \Dij \ in 
meV for Fe, Co and Ni clusters deposited on Ir(lll), Pt(lll) 
and Au(lll). The icons in the left column indicate the corre- 
sponding cluster geometry as well as the cluster sites i and j, 
respectively. The last line gives \Dij \ for nearest neighbouring 
sites within the full monolayer (ml). 



ij 


ll 


Fe 

rt 


All 


ii 


Co 

_r t 


Au 


ll 


Ni 

it 


Au 


mm 


1 1 7 
1.1 1 


ft n9 

o.yo 


1 ft 1 

l.Ol 


9 A Q 
O.40 


£. AT 

0.4 ( 


O Q A 


n oft 


A OA 

U.z4 


A A(\ 

u.4y 




A ftn 

4.0U 


ft i ft 
0.10 


A 7H 
U. (U 


i nft 
1.U0 


1 IT 
1.1 ( 


4.1o 


n 9n 


1 A G 

1.4o 


A AA 

u.yy 





n ciA 
U.94 


ft 9 1 

O.ol 


o a g 
Z.4o 


C 7ft 

0.(0 


o.ol 


Q ftft 
O.OO 


n O/i 
U.z4 


O Q G 
Z.OO 


1 QO 
l.OZ 


omm 


1.83 


k ft a 
0.64 


O T7 


o nft 
2.06 


Q G 1 

o.ol 


1 07 

L.Z{ 


0.40 


a ao 
U.DO 


n ao 
U.4o 


mom 


9 ftO 
O.OZ 


1.19 


Q 7Q 
O. (O 


n qq 
U.oo 


1 T(\ 

i. / y 


O A1 

z.yi 


n /i n 
U.4U 


1 OA 

l.zU 


1 AO 

1.4z 


& 




A Qft 
4.00 


ft no 
O.Uz 


1 CO 

l.OZ 


o on 
Z.Zi) 


A fiT 
U.O ( 


1 79 
1.(6 


n i /i 
U.14 


A Q7 

u.y ( 


1 A G 

1.4o 


.°o 
• 


1 7ft 
1.(0 


C ft A 

0.04 


O Oft 

Z.ZO 


A 99 

4.oo 


A OA 

4.Z4 


z.oo 


n on 
U.zU 


1 AA 

i.uy 


1.0 ( 


& 




O £ A 
Z.04 


r Qn 
O.oU 


A HQ 


i in 
1.19 


i on 

1.03 


O 1 A 

Z.14 


n oft 


A OQ 
U.ZO 


A QA 

U.oy 


A 

• 


1.00 


n ft/i 
U.04 


A 7C 


n 9 a 
U.o4 


1. r O 


A AA 

u.yy 


n nn 

u.uy 


A GQ 

U.oo 


A 1 A 

u.iy 


o 


0.(0 


O.ol 


A ftQ 
U.OO 


O (17 


ft aq 
O.Uo 


T AQ 

( .yo 


n £ /i 
U.04 


A C7 

U.o / 


A Oft 

U.ZO 


..°o 

o 


Z.Zo 


D.LZ 


1 oc 
l.OO 


A CI 

4.0l 


g ftQ 
o.oo 


Q OA 

o.zy 


n £o 
U.Oz 


1 AO 

1.4z 


A 7Q 
U. Id 


oA 

o 


A ft7 
4.0 I 


ft (17 

o.y ( 


o gi 
z.ol 


9/19 

o.4o 


ft OA 

0.z4 


K QO 

O.oz 


n i 
U.01 


A QQ 
U.OO 


A AA 

u.yu 


•o°o 
• 


l.o» 


K A ft 

0.40 


1 /I G 

1.4o 


( .00 


A AC 

y.yo 


4. (1 


n Kn 
U.OU 


l.UD 


1 A Q 

1.4o 





l.oU 


1 HQ 

l.Uo 


O OA 
Z.Z4 


n Q9 
U.Oo 


n qa 
U.oU 


1 AC 

i.yo 


n 97 
U.o ( 


U.oo 


A OQ 
U.ZO 


•A 

O 


O 9/1 


n 77 
U. ( ( 


A Qft 

U.oo 


1 At\ 

1.4U 


z.zl 


Oil 

o.ll 


n i ft 
U.10 


U.Ol 


A Af\ 

U.4U 


.. 

00 


f; fin 
O.OU 


K Oft 

O.zO 


O AA 
Z.UU 


o 9n 
Z.oU 


1 QO 
l.OZ 


G QA 

o.oy 


U.10 


A QQ 
U.OO 


1 OI 

l.zl 


.# 


3.33 


3.64 


1.19 


1.14 


1.47 


1.34 


0.26 


0.54 


0.34 


c# 


2.51 


5.22 


1.43 


3.83 


1.54 


3.81 


0.39 


0.42 


1.19 




2.90 


4.90 


1.21 


1.46 


2.58 


2.40 


0.31 


0.28 


0.76 


.# 


0.79 


0.35 


0.73 


0.40 


1.73 


1.10 


0.17 


0.40 


0.49 




1.49 


0.46 


1.45 


1.12 


1.31 


1.74 


0.35 


0.74 


0.27 




1.00 


4.99 


1.08 


5.01 


1.86 


1.99 


0.45 


0.26 


0.72 




3.72 


6.00 


2.55 


1.44 


2.88 


1.56 


0.26 


0.23 


1.22 




1.28 


0.43 


0.97 


0.57 


1.20 


0.80 


0.10 


0.27 


0.39 


•• 


1.28 


0.43 


0.97 


0.57 


1.20 


0.84 


0.10 


0.29 


0.38 




0.58 


0.43 


1.70 


1.97 


1.56 


0.86 


0.12 


0.33 


0.33 


# 

ml 


3.81 
4.17 


2.52 
2.54 


2.55 
0.66 


1.48 
1.33 


1.75 
1.68 


1.67 
0.71 


0.13 
0.13 


0.35 
0.20 


0.45 
0.04 



collinearity will not be a consequence of the frustration 
between the magnetic and geometric order but rather 
will follow from the influence of spin-orbit effects on the 
exchange coupling, as manifested by the DM interaction. 



D. Magnetic Anisotropy Energy 

The magnetic anisotropy energies (MAE) per atom for 
all investigated clusters are compiled in Table IIVI Posi- 
tive MAE values denote an out-of-plane anisotropy while 
negative MAE values correspond to an in-plane mag- 
netic easy axis. Fe clusters on Pt(lll) and Au(lll) 
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TABLE IV: Magnetic anisotropy energy (MAE) per atom for 
Fe, Co and Ni clusters deposited on Ir(lll), Pt(lll) and 
Au(lll). The icons in the left column indicate the corre- 
sponding cluster geometry. The positive (negative) values of 
the MAE (in meV) correspond to an out-of-plane (in-plane) 
magnetic easy axis. For the full monolayers the total MAE 
AE is decomposed into its dipolar part AEdd and its mag- 
netocrystalline part AE soc . The latter is further decomposed 
into contributions that originate from the monolayer (AE^ C ) 
and the substrate (Ai?£f c ), respectively. For the deposited 
clusters we found A£dd « and AE BOC « AE™ C . 
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show always an out-of-plane MAE whereas all other clus- 
ter substrate systems exhibit a rather nonuniform be- 
haviour of their MAE with varying cluster size or geom- 
etry. This complex behaviour arises from the fact that 
already tiny changes in the electronic structure can cause 
large changes in the MAE. This can be seen again for ex- 
ample in the case of the compact trimers where one can 
observe a dramatic dependence of the MAE on their po- 
sition with respect to the substrate, i.e. depending on 
whether a substrate atom is underneath the cluster cen- 
tre or not. 

All dimers and linear trimers with in-plane MAE have 
their magnetic easy axis fixed along the cluster axis which 
is a result of the strong azimuthal MAE in these systems 
being in the order of 1-4 meV per atom. For compact 
symmetric clusters as well as for the full monolayers there 
remains only a very small azimuthal MAE in the order 
of /LteV being thus negligible. 

When evaluating the MAE by means of the torque 
method, contributions stemming from all individual 
atomic sites of the system are added together. One can 
therefore technically identify which portion of the MAE 
comes from the adsorbed atoms and which portion comes 
from the substrate atoms. We found that the contribu- 



tion coming from the substrate is negligible in the case of 
clusters while it can be significant in the case of complete 
monolayers (e.g. up to 45% of the total value for the Co 
monolayer on the Pt(lll) substrate). This is plausible 
given the fact that for monolayers, the substrate atoms 
are subject to interaction with a larger number of ad- 
sorbed atoms, meaning that their spin-polarization will 
be stronger and more robust than in the case of small 
clusters, contributing thereby more significantly to the 
MAE. At the same time, one has to bear in mind that 
energy is not an extensive quantity and that any decom- 
posion of the MAE into parts has only a limited signifi- 
cance. 

Concerning the dipole-dipole or shape MAE contri- 
bution, for clusters it is negligible while for complete 
monolayers it attains appreciable values of -0.19 meV and 
-0.09 meV per atom for Fe and Co monolayers, respec- 
tively. Moreover, we find that the substrate as well as 
the dipole-dipole contribution to the MAE is negligible 
for clusters, whereas for monolayers both contribution 
are much more important. 



E. Comparison with other works 

As already mentioned in the introduction it is not 
always straightforward to directly compare theoretical 
LSDA results obtained by different computational ab- 
initio implementations due to differences in the trunca- 
tion of the wavefunction or Green's function etc. as well 
as different technical issues and approches as for example 
the implementation of spin-orbit coupling as perturba- 
tion, the use of a supercell vs. embedding techniques or 
approximations in the description of the effective poten- 
tials and so forth. All this can affect the obtained numer- 
ical results, especially for sensitive magnetic quantities 
like for instance orbital magnetic moments and magnetic 
anisotropy energies. 

Among the cluster/substrate systems discussed in this 
article, only Fei and Co„ on Pt(lll) have been stud- 
ied extensively by other groups and we find for these 
systems that our spin magnetic moments agree quantita- 
tively well with the corresponding /i sp i n values given in 
Refsi^i^ir— using identical geometries. The same is true 
for Fe 3 on Pt(lll)^, Fei and Coi on Ir(lll)21 as well as 
for the monolayer systems Fe/Ir(lll) 27 i , Co/Ir(lll)2i, 
Fe/Pt(lll), Co/Pt(lll) 31 and Co/Au(lll) 37 . Regard- 
ing the values of /z or b and the MAE, however, the agree- 
ment is in general less good, i.e. only qualitative or 
worse, for the above mentioned reasons. As already anal- 
ysed by Sipr et al/ 3 ^, methods which rely on a super- 
cell approach 3 ^— produce always significantly higher in- 
duced spin magnetic moments within the substrate atoms 
when compared to methods which apply embedding tech- 
niquesi^— . 
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IV. SUMMARY AND CONCLUSIONS 

The evolution of the spin and orbital magnetic mo- 
ments of the investigated 3c? transition metal clusters on 
5d noble metal surfaces mostly follows common trends 
and patterns that can be understood by considering the 
coordination numbers of atoms in the clusters and the 
polarizability of the substrate. The average /i sp i n values 
decrease nearly monotonously with the number of atoms 
in the cluster being at variance with trends observed for 
free clusters^. Our results show that /i or b may strongly 
depend on the position of the cluster with respect to the 
surface atoms, as demonstrated in particular for the tri- 
angular 3-atom clusters on Pt(lll) and Au(lll). The 
magnetic moments for Ni clusters on Ir are smaller than 
one would expect judging from the trends for the other 
cluster/substrate combinations. Moreover, they depend 
wildly on number of atoms in the cluster and their small- 
ness is compatible with the fact that the peak in minority 
DOS is below E F . 

Apart from Ni„/Ir(lll) all clusters show a strong ferro- 



magnetic isotropic exchange coupling exceeding the cor- 
responding bulk values of standard bcc Fe, hep Co and 
fee Ni. In addition, there are also strong anisotropic 
DM interactions present revealing the intrinsic tendencies 
towards noncollinear magnetism in these systems. Fi- 
nally, the magnetic anisotropy energies can be very large 
for some cluster/substrate or surface/substrate combina- 
tions, but unfortunately, there are no clear trends visible 
that would allow any straightforward anticipation of this 
sensitive quantity. 
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